On the Detectability of the Cosmic Dark Ages: 
21 -cm Lines from Minihalos 
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^ ' Abstract. In the standard Cold Dark Matter (CDM) theory of structure formation, viriahzed mini- 

halos (with Tvir < lOjOOOK) form in abundance at high redshift (z 6), during the cosmic "dark 
' ages." The hydrogen in these minihalos, the first nonlinear baryonic structures to form in the uni- 

verse, is mostly neutral and sufficiently hot and dense to emit strongly at the 21 -cm line. We cal- 
culate the emission from individual minihalos and the radiation background contributed by their 
combined effect. Minihalos create a "21 -cm forest" of emission lines. We predict that the angular 
^ I fluctuations in this 2 1 -cm background should be detectable with the planned LOFAR and SKA radio 

arrays, thus providing a direct probe of structure formation during the "dark ages." Such a detection 
will serve to confirm the basic CDM paradigm while constraining the background cosmology pa- 
rameters, the shape of the power-spectrum of primordial density fluctuations, the onset and duration 
^— ^ ■ of the reionization epoch, and the conditions which led to the first stars and quasars. We present 

I results here for the currently-favored, flat ACDM model, for different tilts of the primordial power 

. spectrum. These minihalos will also cause a "21 -cm forest" of absorption lines, as well, in the spec- 

' trum of radio continuum sources at high redshift, if the latter came into existence before the end of 

'q^I reionization. 

o 

^ : INTRODUCTION 

No direct observation of the universe during the period between the recombination epoch 
^ . at redshift z — 10^ and the reionization epoch at z ^ 6 as yet been reported. While a 

number of suggestions for the future detection of the reionization epoch, itself, have 
been made, this period prior to the formation of the first stars and quasars - the cosmic 
"dark ages" - has been more elusive. Standard Big-Bang cosmology in the CDM model 
predicts that nonlinear baryonic structure first emerges during this period, with virialized 
halos of dark and baryonic matter which span a range of masses from less than IO'^Mq to 
about 10^ Mq which are filled with neutral hydrogen atoms. The atomic density and 
kinetic temperature of the gas are high enough that collisions populate the hyperfine 
levels of the ground state of these atoms in a ratio close to that of their statistical 
weight (3:1), with a spin temperature Ts that greatly exceeds the excitation temperature 
= 0.068 IK. Since Ts > Tqmb, the temperature of the Cosmic Microwave Background 
(CMB), as well, for the majority of halos, these "minihalos" can be a detectable source 
of redshifted 21 -cm emission. The direct detection of minihalos at such high redshift 
would be an unprecedented measure of the density fluctuations in the baryons and of the 
total matter power spectrum at small scales. The results presented here are described in 



more detail in Jill 0]. 



21-CM EMISSION AND ABSORPTION 

The 21 -cm emission from a single halo depends upon its internal atomic density, tem- 
perature, and velocity structure. We model each CDM minihalo here as a nonsingular, 
truncated isothermal sphere ("TIS") of dark matter and baryons in virial and hydro- 
static equilibrium, in good agreement with the results of gas and N-body simulations 
from realistic initial conditions LIl,|2|,|7i]. The minihalos which contribute significantly to 
the 21-cm emission span a mass range from M^i^ to Mmax which varies with redshift. 
A/min is close to the Jeans mass Mj of the uncoUapsed IGM prior to reionization, while 
Mmax is the mass for which Tyir = lO^K according to the TIS model Typically, at 
z = 9, Mmin ~ few X 10^ Mq and Mmax ~ few x IO^Mq, depending upon the cosmological 
model. 

Our results for individual minihalos are summarized in Figure [1^. Line profiles of 
different minihalos along the same line of sight should not typically overlap. The proper 
mean free path for photons to encounter minihalos in ACDM is 160 kpc at z = 9 yD, 
corresponding to a frequency separation AVsep ~ 0.1 MHz ^ AVeff ^ lOkHz. These 
results predict a "21-cm forest" of minihalo emission lines. At z = 9,_for example, 
there are about 160 lines per unit redshift along a typical line of sight|5]. Detecting 
the stronger lines would require sub-arcsecond spatial resolution, ~ IkHz frequency 
resolution, and ~nJy sensitivity. SKA is expected to have sufficient resolution for such 




FIGURE 1. (a) (left) 21-cm emission line profiles for individual minihalos of different mass at z = 9. 
Differential antenna temperature 57], [K] versus emitted frequency Vem, from detailed radiative transfer 
calculation (lower curves) and optically thin approximation (upper curves). Optical depth is important, 
particularly for smaller halos. (b) (right) 21-cm absorption line profiles for individual minihalos of 
different mass at z = 9. Transmission factor versus received frequency Vrec at z = for line of sight 
through minihalo center (i.e. zero impact parameter). 



observations, but probably not sufficient sensitivity. 

The minihalos which produce the "21 -cm forest" of emission lines described above 
have appreciable optical depth to 21 -cm absorption. The same minihalos will therefore 
produce a "21 -cm forest" of absorption lines, too, in the spectra of radio continuum 
sources if the latter are discovered at high redshift (z > 6), prior to the end of the 
reionization epoch. We plot illustrative 21 -cm absorption line profiles for absorbing 
minihalos of different masses at z^bs = 9 in Figure[lj5. These absorption features should 
be observable provided sufficiently bright (~ 1 mJy) background source exists. 



21-CM RADIATION BACKGROUND 



The beam- averaged differential temperature dT}, within a given beam of angular size 
AQbeam is Calculated by integrating our results for individual minihalos over the Press- 
Schechter mass function of minihalos sampled by the beam. We consider the currently- 
favored fiat ACDM model (Qq = 0.3, = 0.7, COBE-normalized, Qbh^ = 0.02, h = 
0.7), with primordial power spectrum index Hp = 0.9, 1, and 1.1. 

In principle, the variation of 5Ti, with observed frequency should permit a discrimina- 
tion between the 21 -cm emission from minihalos and the CMB and other backgrounds, 
due to their very different frequency dependences. However, the average differential 



B 



- 1 \ 

r 1.1 

" To ■ — 




\? 




579 
ACDM 




"-^^^^\^ 




_ ^u= \ MHz 








- z = 7 


•.a 


'■o 




LOJI 1 1 ll| 






M 1 M 


— -CClZ — 




>I<1 




0.9 




^ -.CI \'ci 












[ z = 8.5 


■ o 






" 1 









W X) 

^ 1 




2 

g 1.5 



•o 0.5 



r-L^A ' } 


.-■■'acdm - 


r 7^ — ~— ii 


\ MHz- 






-r^xm^.s ><;^ 




y \^ ... ■ iOOQh.S 




beam 




1 / 1 1 1 1 ; 1 1 1 
/lOO" 

: — — 1.1 : 


\. . ^\ 1.0 




\^ . iOOQi' \^ 




beam | ^ 





1 10 100 10 15 20 

A6(arcmin) z 

FIGURE 2. (a) (left) Predicted 3-(7 differential antenna temperature fluctuations at z = 7 (Vrec = 
177.5MHz, top panel) and z = 8.5 (Vrec = 150MHz, bottom panel) for bandwidth AVobs = IMHz vs. 
angular scale A0beam for ACDM models with tilt rip — 0.9, 1.0, and 1.1, as labeled (solid curves). Also 
indicated is the predicted sensitivity of LOFAR and SKA for a confidence level of 5 times the noise 
level after integration times of 100 hr (dashed lines) and 1000 hr (dotted lines), with compact subaperture 
(horizontal lines) and extended configuration needed to achieve higher resolution (diagonal lines), (b) 
(right) Predicted 3-(7 differential antenna temperature fluctuations at A0beam = 9' (top panel) and 25' 
(bottom panel). Symbols have same meaning as in the left panels. Letters "L" and "S" correspond to 
LOFAR and SKA, respectively. 



brightness temperature of the minihalo background is very low and its evolution is fairly 
smooth, so such measurement may be difficult in practice with currently planned in- 
struments like LOFAR and SKA. The angular fluctuations in this emission, on the other 
hand, should be much easier to detect. The amplitude of q-o angular fluctuations (i.e. q 
times the rms value) in the differential antenna temperature is given in the linear regime 
by {dT^y/'^ / 5Th = qb{z)Op, where Op is the rms density fluctuation at redshift z in a 
randomly placed cylinder which corresponds to the observational volume defined by the 
detector angular beam size, A^beam^ and frequency bandwidth, AVobs^ and b{z) is the bias 
factor which accounts for the clustering of rare density peaks relative to the mass. 

Illustrative results are plotted for 3-C7 fluctuations as a function of A0beam for z = 7 
and 8.5, in Figure EK, along with the expected sensitivity limits for the planned LOFAR 
(300 m filled aperture) and SKA (1 km filled aperture) arrays. We plot in Figure |2h the 
predicted spectral variation of these fluctuations vs. redshift z for illustrative beam sizes 
of A0b earn — 9' and 25'. These 3-c7 fluctuations should be observable with both LOFAR 
and SKA with integration times of between 100 and 1000 hours. For a 25' beam, for 
example, 3-a fluctuations can be detected for untilted ACDM by both with a 100 hours 
integration for z ~ 6 — 7.5 and a 1000 hours integration for z ^ 11.5, while for a 9' beam, 
SKA can detect them after 100 hours for z ^ 9 and after 1000 hours for z ^ 13. 
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